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Abstract 

We study the real time correlators of scalar glueball operators for Yang-Mills the- 
ory at finite temperature in flat space. The analytic structure of the frequency space 
propagator in perturbative field theory is seen to be qualitatively different to the strong 
coupling results that may be obtained from perturbations about AdS black hole space- 
times: we find branch cuts rather than poles. This difference appears to persist away 
from the strict zero and infinite coupling limits, possibly suggesting a phase transition 
in large N thermal J\f = 4 SYM theory as a function of the 't Hooft coupling. 
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1 Introduction and summary 

The AdS/CFT correspondence relates gravitational and field theories in complementary 
regimes of tractability [1]. This fact is both an important success and a significant limitation 
of the correspondence. On the one hand, such complementarity has allowed insight into 
the dynamics of strongly coupled gauge theories. However, if one wishes to study quantum 
gravity processes involving large curvatures in the interior but with asymptotically weak 
curvatures, such as evaporation of black holes, then the appropriate dual description is a 
strongly coupled field theory. Thus neither classical bulk physics nor perturbative boundary 
physics is able to provide significant insight into these processes. 

It is therefore an important recent observation that gravitational physics in asymp- 
totically AdS space with large curvatures everywhere appears to posses a similar phase 
structure to gravity with asymptotically weak curvatures [2, 3, 4, 5, 6]. In particular, there 
is a Hawking-Page phase transition [7] and a description in terms of an effective potential 
with two minima and a local maximum that may be interpreted as thermal AdS space, the 
large AdS black hole and the small AdS black hole, respectively. The emerging picture sug- 
gested that perturbative field theory could provide an inroad for studying quantum gravity 
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processes that were at least qualitatively similar to those of real physical interest. Thus for 
instance one might hope that generalising the work of [3, 4, 5] to consider out of equilibirum 
real time thermal field theory on S" 3 x R could provide a framework for studying black hole 
evaporation. 

The main thrust of this work is to exhibit in some detail one aspect at least in which 
black holes in strongly curved AdS space are qualitatively different to the more familiar 
weakly curved black holes. One consequence of these differences will be that the strongly 
curved black holes are not particularly black. Fields decay exponentially in weakly curved 
AdS black hole backgrounds [8, 9], with a timescale set by the black hole temperature. In 
free field theory on flat space, we see that the decay is generically only power law whilst 
including weak interactions results in a power law times a slow exponential decay over the 
plasmon damping timescale l/(ATln(A" 1 )) , where A <C 1 is the 't Hooft coupling. This 
behaviour is in line with recent observations in [10], who noted that black holes dual to 
weakly coupled field theories are not totally absorbent. 

We will study the finite temperature, frequency space retarded propagators of the scalar 
glueball operators TvF^F^ and TvF^F^ at both strong and weak 't Hooft coupling in 
Yang-Mills theory in Minkowski spacetime. The M = 4 supersymmetric Yang-Mills (SYM) 
theory at finite temperature and in flat space is realised in a deconfined phase. In this 
context the operators we study are dual to the dilaton and axion (the RR scalar) fields, 
respectively, in the bulk AdS black hole background. At weak coupling, a gauge theory 
computation indicates that these propagators generically contain branch cuts. At strong 
coupling, the same propagators have poles and no branch cuts as found in, for example, 
[11, 12]. We build on previous computations of the strong coupling behaviour to argue 
that these poles persist under a 1 corrections without the introduction of branch cuts. This 
difference in analytic structure is directly related through a Fourier transform to the differing 
time decays for the fields mentioned in the previous paragraph. Interestingly, we find that 
the weak coupling frequency space propagators bear striking similarities to corresponding 
objects computed in topological black hole backgrounds such as the BTZ black hole. 

The structure of the paper is as follows. In section 2 we compute precisely the free 
Yang- Mills glueball propagator on R 1 ' 3 , giving both the time dependence and frequency 
space results. We also consider the effects of perturbative corrections to these results, 
noting in particular the resummations that are necessary due to hard thermal loops and 
collisions [13, 14, 15]. In section 3 we recall the prescription for computing finite temperature 
correlators from bulk black hole backgrounds [11, 16]. We use this prescription to compute 
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and plot the retarded propagators in frequency space. In section 4 we prove that at least 
some of the strong coupling poles persist under a! corrections whilst branch cuts do not 
appear in their vicinity. Finally, section 5 contains a discussion in which we speculate on 
possible interpretations of our results, such as the existence of a phase transition in large 
iV thermal N = 4 SYM theory at an intermediate 't Hooft coupling, and consequent future 
directions for research. Appendix A explains the effect on the correlators of including the 
thermal masses and widths, and Appendix B contains computations of free thermal glueball 
correlators on 5 3 x R. 

2 Thermal glueball propagators on R 1 3 

We begin by considering SU (N) J\f = 4 supersymmetric Yang-Mills (SYM) theory at finite 
temperature T in four dimensional Minkowski spacetime and at small 't Hooft coupling, 

A = g 2 YM N « 1. (1) 

In this section we will analyse the finite temperature, real time behaviour of two point 
correlators of the marginal, 50(6) singlet scalar glueball operators 

G = TrF^F^ , g = TrF^F^. (2) 

These are chiral operators in the = 4 theory. In the Euclidean or static description 
wherein the thermal field theory is formulated on R 3 x 5 1 these operators couple to scalar 
glueball states of the three dimensional effective theory. 

The temporal decay of fluctuations or correlations in a weakly coupled Yang- Mills plasma 
at finite temperature will generally occur at several distinct time scales. At short time scales 
set by the inverse temperature, t ~ 1/T, the behaviour is completely determined by free 
field theory. Intermediate and late time behaviours are controlled by other physical scales 
such as the thermal gluon mass scale (\/AT) _1 , the gluon damping rate (AT) -1 , etc. [14]. 

2.1 Short times or the free theory 

For short time scales t ~ T" 1 , and weak coupling, thermal correlators can be approximated 
by straightforward diagrammatic perturbation theory. Thus, at leading order we only need 
the quadratic terms in the Yang- Mills action and the glueball operator. Furthermore, we can 
focus attention solely on the transverse, physical degrees of freedom, since the longitudinal 
mode only mediates an instantaneous Coulomb interaction and does not affect the dynamics 
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at any finite time t > 0. We also therefore don't need to worry about ghost fields. The time 
ordered two point function for the transverse gauge field is 

<i4?(t,p)A}(o, P )> = D>/ ab \t) e(t) + D< ab] \t) e(-t) (3) 

where p is the spatial momentum and [17, 18] 

D> iab) (t) = D<<*> (-t) = i (dy - g [(1 + n(£0)e"^ + n(£)e^] . (4) 

Here E = p = |p| and n(-E) is the Bose-Einstein distribution 

n ( E ) = jnnT[> /3=|- (5) 

The glueball two point function with external momentum p is then straightforwardly worked 
out using Wick contractions to be 



(G p (t)g p (o)) = (G p (t)G p (o)) = J d 3 x e - i p- x (g(i,x)g(o)) 
= -£ / IpIIp - p1 [(p/2(p ~ pr + [p/ - (p ~ p0]2) F(|p/| ' t)F(|p " p/| ' t} 



+2p / -(p-p / ) J F(|p / |,t)F(|p- P , |,t)J , (6) 

where 

F(E, t) = (n{E) + l)e-' lEt + n(E)e iEt . (7) 

The Q and Q glueballs have the same free theory thermal two point functions. We will write 
Q from now on for notational convenience. 

We first consider the homogeneous mode with no external momentum, p = 0. The 
relaxation of long- wavelength fluctuations or 'soft' modes withp <C T~ x can be non-trivially 
altered at late times by interactions, as we discuss below. Nonetheless in the limit of zero 
coupling, or for suitably short time scales in the interacting theory, the correlator at zero 
momentum, with p = in (6), becomes 

o at2 roo 2x „— i2xr i J!2xt 

(OW(0)) = -^j o d XX * e e [ex _ + x) l , (8) 

where we have introduced the dimensionless time r = tT and the variable of integration 
x = \p'\/T. The integral may be evaluated to give 
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(0°(r)g°(0)) = ^T 5 



8 ^ + (2T + ^ 



coth(2vrr) . (9) 



As r — > there is a 1/r 5 divergence, which is the zero temperature power law behaviour. 
Importantly, this expression exhibits an exponential decay on time scales r ~ 1 or t ~ T _1 . 
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Using the above result we obtain the retarded propagator which is of particular physical 
interest and is given by 

N 2 rl 4 

(g°(t)g°(o)) R = -i&(t)([g (t),g°(o)}) = —t @(t)^ coth^rt) . (10) 

The retarded propagators naturally occur in the framework of linear response theory wherein 
one probes the plasma with a weak disturbance and studies the response of the medium in 
a linearised approximation. 

In the AdS/CFT correspondence, retarded thermal propagators of the boundary field 
theory are known to encode physical properties of AdSs black hole backgrounds, at least in 
the classical gravity approximation which is dual to strongly coupled field theory. Specif- 
ically, poles of retarded propagators in the frequency domain correspond to quasinormal 
frequencies of AdS black holes [19, 11]. The scalar glueball operators g and g are dual to 
the dilaton and axion in the bulk, respectively. The associated retarded correlators at strong 
coupling encode quasinormal frequencies of massless minimally coupled scalar perturbations 
about the AdSs black hole background. We will return to the gravitational computation of 
correlators in later sections of this paper. 

A characteristic feature of scalar perturbations about AdS black hole spacetimes is 
that they always decay exponentially [8, 9]. The exponential decay (10) that we have 
encountered in the free theory might seem to be consistent with this feature, although it 
is in a regime where the dual geometry is very strongly curved. However, the exponential 
decay of quasinormal perturbations in AdS black hole spacetimes represents the decay 
of perturbations in the boundary thermal theory and its consequent approach towards 
equilibrium. The exponential decay that we see for the homogeneous mode in the free 
theory is not, strictly speaking, an approach to equilibrium since scattering and collisional 
processes are absent in this limit. The decaying transient that we see here arises from 
interference effects between the different modes in (8). In fact we will see shortly that all 
the inhomogeneous modes of (6) in the free theory have oscillating power law decays instead, 
a fact which will have a simple physical explanantion. 

The qualitative similarities and dissimilarities between the field theory at weak and 
strong 't Hooft couplings are naturally captured by the analytic structure of the retarded 
propagator in frequency space. Fourier transforming equation (10) we find 



-ioj\ 1 /2vrT\ 2 2 /2vrT\ 4 2ttT 



UJ 



(11) 



where tp(z) = T'(z)/T(z), the logarithmic derivative of the Gamma function. In the zero 
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temperature limit, using Stirling's approximation for T(z) at large z, we have 



i-K sgn u) + ■ ■ ■ as — — > 



oo , 



(12) 



which is the known zero temperature result (see e.<?.[ll]). 

The Fourier transformed homogeneous correlator (11) is analytic in the upper half plane 
3m u> > 0. This is always true for retarded propagators with our conventions. On the other 
hand the function has an infinite set of simple poles along the lower imaginary axis at 



These poles, illustrated in figure 1, lead to the exponential decay in time of (10). This 
analytic structure at zero 't Hooft coupling will of course be modified by a small nonzero 
value of A and we will discuss the modifications in more detail below. Perhaps unsurprisingly, 
the analytic structure at A = bears little resemblance to that at A — > oo wherein the 
retarded correlator has an infinite number of simple poles at locations that at large n tend 
to uj^ ± 2irTn(l =p i) where lo^ ~ ±1.2139 — 0.7775i [11]. The strong coupling poles 
are sketched in figure 2. The strong- weak coupling discrepancy will be substantially more 
severe for the inhomogeneous modes. 

It is interesting to note that the frequency space propagator (11) is, up to a dimension 
dependent phase space factor, precisely the function encountered for the propagator at zero 
spatial momentum in the general BTZ black hole background dual to a two dimensional 
CFT [11]. Drawing further on this analogy, perhaps this suggests that like the BTZ black 
hole, the strongly curved five dimensional background is a kind of topological black hole 
(c.f. higher dimensional AdS space with global identifications in for instance [20, 21, 22]). 
From the free field theory perspective, the similarity to a two dimensional CFT and the 
BTZ black hole arises simply because free field correlators are completely determined by 
the scaling dimensions of the operators in question and the effect of finite temperature is to 
require an infinite sum over images that renders the correlators periodic in imaginary time. 

The inhomogeneous contribution to the propagator (6), with spatial momentum p, may 
also be calculated. Introducing the dimensionless momentum k = p/T, the retarded prop- 
agator is given in terms of the integral 



= —4mTn 



n= 1,2... 



(13) 



<S p (t)^(0)) r = 




( e (u+v)/2 _ l^ e (u-v)/2 _ 



(u 2 -k 2 ) 2 (l-e u ) 



sin(wr) 




sin(vr) 



(14) 



( e (u+«)/2 - i)( e (u-«)/2 _ !) 
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Figure 1: Poles and branch cuts of the free theory propagator G p r (uj). Top left is the full 
result with momentum p. Top right is the homogeneous mode with p = 0. Bottom left is 
the zero temperature limit with momentum p. 



The double integral (14) may be evaluated, partly using contour integration, to give 



(Q p (t)G p (0))R 



sin pt 



7T coth(2vrtT) 



sin pt — pt cos pt 



(15) 



r p v dt 2 J [ t v y 2t 2 T 

As £ — ► oo this propagator decays as an oscillatory power law rather than exponentially. 
The physical reason for this is that these modes carry a conserved momentum which in the 
free theory cannot dissipate into other modes. In the strongly coupled theory (A — ► oo), 
dual to the classical black hole, the momentum can quickly dissipate into other modes and 
so one expects an exponential decay in this indeed occurs. A small but nonzero A 

will also modify this power law decay at intermediate and late times as we discuss below. 
The propagator (15) may be Fourier transformed to frequency space to yield 

.irT • ' ' ' 1 

~2p~ 



N 2 

-2tt 2 T 2 



+- 



7T Z 
N 2 



1 

2 + 
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7T^ 
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16vr 4 T 4 
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-i(u>+p) \ 
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( -i(oJ-p) \ 
\ 4wT J 



(16) 



3 ' ' 15 

It is easily seen that (16) reduces to (11) in the vanishing momentum p — > limit. In the 
zero temperature limit we also recover the wellknown Lorentz invariant result, providing a 
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quick check of (16) 

^M = ^(p 2 -^ 2 ) 2 log^. (17) 

The central point to note here is that the frequency space result (16) has only branch cuts 
and no poles, as illustrated in figure 1. The infinite set of simple poles of the zero momentum 
correlator (11) have smeared out into branch cuts stretching between uj = —4niTn — p and 

uj = —AiriTn + p with n = 1, 2, In the time domain, these branch cuts centred around 

points on the negative imaginary axis translate into exponential falloff accompanied by 
power law and oscillatory corrections. In addition to these, there is also a branch cut on 
the real axis between u> = — p and oj = +p and it is this branch cut which is responsible for 
the leading oscillatory power law behaviour in time in (15). 

These results may be contrasted with the analytic structure emerging from dual gravi- 
tational computations that describe the strongly coupled regime, e.g. [12, 11]. The strong 
coupling results are sketched in figure 2, displaying an infinite set of simple poles correspond- 
ing to the quasinormal frequencies of the AdSs black hole background. These asymptote 
to u)q{p) ± 2irTn{\ =p i). Note that both the real and imaginary parts of the quasinormal 
frequencies grow linearly with n = 1, 2, 
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Figure 2: Poles in the strongly coupled propagator [12]. 



Thus we can see that correlators in the A = theory exhibit features that are qualita- 
tively different to those of the A — > oo theory. As we will argue in the following section, 
moving slightly away from the A = point does not change this qualitative picture dras- 
tically, although it does affect the late time behaviour of the correlator (15) causing it to 
decay exponentially. While branch cuts in frequency space induce power laws in the time 
domain, the endpoints of the branch cuts dictate the exponential decay constants and/or 
oscillatory behaviour. For fixed p the real parts of the branch point locations of (16) are 
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fixed while the imaginary part increases with n. Following our previous comments concern- 
ing parallels with the BTZ black hole, it is tantalising to note that quasinormal frequencies 
of topological AdS black holes appear to have a similar property, namely that for fixed 
spatial momentum their real parts are constant while the imaginary parts take an infinite 
set of discrete values [23] . 

Finally, we have already noted that the operator Q (x, t) in (2) has exactly the same 
two point function as Q in free field theory. This operator has the quantum numbers of 
the h scalar glueball, and at strong coupling it is dual to the axion (the RR scalar) in 
supergravity. We will see shortly that this implies that also at strong 't Hooft coupling 
the operators Q and Q share the same correlators. At finite values of A however, we expect 
them to differ quantitatively. 

2.2 Late times and weak coupling A ^ 

We now turn to the effects of nonzero 't Hooft coupling A on the late time behaviour of 
the correlators above. In the interacting N = 4 SYM theory, the correlators of interest 
receive higher order contributions from gluon self-interactions and from scalar and fermion 
matter fields in the adjoint representation of the gauge group. Since the thermal bath breaks 
SUSY and conformal invariance, in what follows it is sufficient to focus on the known general 
qualitative aspects of time dependent correlators in a weakly interacting Yang-Mills plasma 
in the deconfined phase. In the weakly coupled plasma there is a hierarchy of energy scales, 
(T, \/XT, XT, . . .) corresponding to distinct physical phenomena that dictate the temporal 
behaviour of correlators over the associated characteristic time scales. 

In order to go beyond the leading perturbative result it is necessary to implement the 
finite temperature improvement of resumming the hard thermal loops [13]. Essentially this 
means that for the internal lines with momenta p' or p— p' in the above one loop calculation, 
we should include the thermal self-energy corrections. Specifically, for the hard momenta of 
order T which dominate the graph, this self-energy correction leads to an effective thermal 
mass vriT ~ \[XT. In addition there is also a higher order effect leading to a thermal 
width, or more precisely the plasmon damping rate for hard gluons, T ~ ATln(A _1 ). This 
resummation has a significant effect on homogeneous correlators or those with 'soft' external 
momentum p < my ~ y/XT. For 'hard' modes, with external momentum p > T, the effect of 
the hard thermal loop resummation remains perturbatively small and does not significantly 
alter the early time oscillatory power law behaviour found in (15). 

First consider the effect of incorporating thermal masses for internal lines. This will 
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shift the threshold branch points at u = ±p to lo = ±yp 2 + 4my. The thermal mass thus 
changes the analytic structure of the homogeneous correlator with p = (11) by introducing 
a branch cut between lo = +2my and lo = —Itxit which leads to oscillatory power law decays 
in time. This change in the analytic structure can be inferred from (6) by setting p = and 
using the Bose-Einstein distribution for a massive particle. We perform this computation 
in Appendix A. 

It is easy to see, following the example in Appendix B of [14], that this inclusion of 
the thermal mass leads to a power law decay with oscillatory behaviour for time scales 
t > (\/AT) _1 , where the oscillations are controlled by the thermal mass tut- Hence the 
short time exponential decay of the homogeneous correlator (10) turns into an oscillatory 
power law on time scales of order (VXT)^ 1 . This power law behavior kicks in when the free 
field correlator has decayed to V\ times its value at t = 0. 

However, including only thermal masses is, strictly speaking, a collisionless approxima- 
tion from the point of view of quasiparticles in the plasma although it includes quantum 
field theory interactions. What happens eventually to all correlators at late times is actu- 
ally determined by scattering and collisions of quasiparticle states. A detailed quantitative 
treatment of this question is beyond the scope of this work. However, much can be learned 
by following the arguments in [14]. Incorporating the thermal widths of intermediate states, 
the imaginary part T of the self energy, makes internal propagators decay on time scales 
of order (ATln(A _1 )) _1 and hence the correlators themselves decay exponentially on these 
time scales. This is a somewhat subtle point as discussed in [14]. The effect of incorporating 
the width, or plasmon damping rate, can in certain cases actually cancel out between differ- 
ent diagrams in resummed perturbation theory so that the exponential decay time scale is 
longer and is dictated by different physics. In particular, this occurs for correlators involv- 
ing conserved currents wherein the time scale for exponential decay is controlled by that of 
large angle scattering ~ (A 2 Tln(A _1 )) _1 and can be seen using an effective kinetic theory 
description [15] 1 . Since the scalar glueball operators are not associated to any conserved 
currents or charges, we expect them to decay on time scales set by plasmon damping which 
is also referred to in [14] as the time scale for 'any-angle scattering'. 

The picture that thus emerges when finite coupling A is considered is that the analytic 
structure of the glueball two point function in frequency space is qualitatively similar to 
1 Correlators for conserved currents in a finite N theory eventually decay via hydrodynamic power law 
tails which are not visible in the large N limit [24]. These 1/N suppressed power law tails correspond to 
non-analyticities in the frequency space correlators at small to. 
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figure 1, with an extra branch cut about the origin along the real axis in the homogeneous 
case. But this description cannot account for the late time exponential decay argued for 
above since the branch cut along the real axis in the frequency plane seen in the lowest 
order approximation would lead to oscillatory power laws at late times. This discrepancy 
suggests the following possibility: 

(1) When all relevant perturbative corrections are accounted for, the branch points on 
the real frequency axis are shifted downwards into the lower half plane by an amount 
~ zATm(A _1 ). A naive calculation illustrating this effect is given in Appendix A. 

This analytic structure for frequency space propagators would be consistent with power 
law tails at intermediate times accompanied by exponential decay times power laws at 
late times. A confimation of this possibility would require a detailed consideration of the 
correlators with full self-energy resummations for the internal lines. Thus we cannot rule 
out the following further possibility: 

(2) The branch cut seen in the lowest order approximation from to = ±2my for the 
homogeneous correlator, and from uj ~ ±p for the inhomogeneous cases, might split 
into a series of poles with perturbatively small separations, which then move off the 
real axis upon inclusion of all perturbative corrections in the large N limit of the 
J\f = 4 theory. 

Whilst we have no evidence in favour of this second possibility, it would reproduce the 
expected intermediate and late time behaviour in the weakly coupled plasma, and would 
also allow a smooth interpolation to the strongly coupled gravity dual picture discussed in 
the following section. 

Regarding the second possibility above, it is not clear how poles might emerge from cuts 
at weak coupling. Indeed, perturbative results for real time response functions in weakly 
coupled plasmas in the deconfined phase generically give rise to branch cuts in correlators 
from zero frequency (see e.g. [25, 26]). 

Our discussion in this section has consequences for bulk physics. The strong coupling 
exponential decay on timescales t ~ T -1 is naturally interpreted gravitationally as being 
due to matter falling into the black hole. The much slower decays we have just described 
at weak coupling therefore suggest that black holes in highly curved AdS backgrounds are 
less absorbent. This phenomenon is consistent with the recent comments in [10]. 
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3 Strong coupling correlators from gravity 

A prescription for computing thermal two point correlators in strongly coupled N = 4 SYM 
theory was developed in [11, 16]. This method has allowed the computation of poles in the 
Minkowski space retarded thermal propagator via bulk quasinormal modes [12, 27], with 
particular success in the hydrodynamic limit [28]. 

In this section we review the bulk prescription and use it to produce plots of the thermal 
propagator as a function of frequency u>. These weak curvature results describe the strongly 
coupled SYM theory. We will compare with plots for the free theory correlators following 
from our results in the previous section, highlighting various qualitative differences. In 
the following section we will show that the strong coupling correlators do not have branch 
cuts, as was implicitly anticipated in [12, 27] and thus confirm figure 2. Our results on the 
absence of branch cuts will extend to include a' corrections. 

For notational clarity when discussing gravitational computations, we shall absorb T 
into the frequency: cj/T — ► to. 

As we wish to compare with field theory on R 1,3 , the bulk is given by the near horizon 
geometry of a stack of non-extremal D3 branes 



ds 2 = ^7T 



-f(z)dt 2 + dx 2 + 



dz 



2 1 



R 2 dnj , (18) 



/(*)J 

where as usual d£l\ is the round metric on a five sphere, R sets the AdS scale and 

/(*) = !- 4' ( 19 ) 

Z H 

with the horizon zh = ft /it. 

The prescription may now be stated as [11, 16] 

Gr{lo) = K^g zz Mz)d z <P*M)\ ZB ^ , (20) 

where K is a normalisation constant. The normalisation is known from the bulk action 
[11] or alternatively can be fixed by the normalisation of the zero temperature behaviour 
u; 4 logu; as u> — > oo. The zero temperature result will agree with the boundary free field 
theory computation because the TrF 2 operator is BPS. The bulk modes (f>(z) are evaluated 
near the boundary zb <C R where they are required to satisfy 4>(zb) = 1- They further 
satisfy incoming wave boundary conditions at the horizon 

M*) ~ (*h ~ z)~ 1 ^ as z^z H . (21) 
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As usual in computing boundary correlators using the AdS/CFT correspodnence, the ex- 
pression (20) will contain volume divergences that correspond to UV divergences in the dual 
theory. In order to meaningfully compare bulk and boundary results, we need to fix a reg- 
ularisation scheme. We do this as follows: Firstly in (20) we minimally subtract off all the 
divergent terms as zb — ► 0. This gives a finite propagator whose only scheme dependence 
is due to the ambiguous subtration of a logarithmic divergence ~ uj 4 log zb / /J> • In order to 
extract scheme independent information we then subtract off from the propagator the terms 
that diverge like uj 4 In oj and uj 4 as uj — > oo. This subtraction corresponds to removing the 
zero temperature result. We will perform the same subtractions on our propagators from 
the previous section. 

The bulk field dual to the glueball operator TrF 2 is the dilaton, whilst the field dual to 
TrFF is the axion. Both of these fields are minimally coupled scalars and therefore have a 
mode equation in the black brane background (18) given by [11, 12] 

,// 1 + u 2 „ w 2 /Att 2 , p 2 /4tt 2 , n 

-4>' + r , ™ <t>--7r — 9^ = °> ( 22 ) 



u(l — u 2 ) u(l — u 2 ) 2 u(l — u 2 ) 

where the coordinate u = z 2 j z\. By solving this equation numerically and varying ui and 
q, we can plot the propagator (20). The real part of the propagator has a more interesting 
structure than the imaginary part, so we plot this. Furthermore, the real part of Gr{uj) 
is symmetric on the real axis, so we restrict plots to uj > 0. It is curious that these two 
operators have exactly the same thermal propagator both at zero and at strong coupling. 

Figures 3 and 4 show the propagator for the homogeneous mode as computed in free 
field theory (11) and from the bulk using (20), respectively. As discussed above, the zero 
temperature result has been subtracted to ensure scheme independence 

AG p r (uj) = G p r {uj) - #(p 2 - uj 2 ) 2 log (p 2 - uj 2 ) . (23) 

Although the plots show some similarity, and curiously both lack an uj 2 growth as uj — > oo, 
they remain fairly distinct, consistent with their differing analytic structure off the real axis. 

Figures 5 and 6 show the effect of including a momentum p in the propagator in free 
field theory and in the bulk, respetively. In this case we have plotted for both positive and 
negative uj to emphasise the features of the plots. 

The main features in the free theory plot of figure 5 are a positive p 4 growth with 
momentum along the uj = axis which is apparent in (16) and a negative uj 2 growth 
which becomes more significant at larger p. The transition between these behaviours occurs 
around uj ~ ±p where the branch points are located. 
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In the bulk computation shown in figure 6 we see that the peak and trough present 
already in the homogeneous case, figure 4, become more pronounced once p is included 
and move linearly like u> ~ ±p as p increases. This is consistent with these peaks being 
resonances of the nearest quasinormal pole to the real w-axis, see figure 2. It was shown in 
[12] that as p is increased the lowest poles move closer to the real axis and tend towards 
ijj = ±p. 

One can numerically check that there is no branch point at, for instance, uj = ±p in the 
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2 

OJ 

27 

Figure 5: The real part of AG p r (oj) as computed in free field theory. 




2tt 

Figure 6: The real part of AG p r (oj) as computed from the bulk supergravity. 

bulk results by computing the propagator along the circles 

oj = ±p + e ie , 9 £ [0, 2vr] , (24) 

and verifying that both real and imaginary parts of the propagator are continuous. In the 
following section we will prove the absence of branch points in some regions of the complex 
oj plane. This argument will then extend to include a' corrections. 
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4 Persistence of poles under a' corrections 

To move away from the strict strong coupling limit, we need to consider string theory a 1 
corrections to the gravity background and to the equations of motion for the fluctuations. 
We will be able to show rigorously that at least some of the poles in the retarded glueball 
propagator that are found in the strong coupling limit persist under a' corrections without 
becoming branch points. The argument does not rely on the precise form of the a' correc- 
tions, which is fortunate as even the leading order corrections to the type IIB string theory 
action are not known completely and hence precise calculations using such corrections are 
not possible in our context, despite some assumptions to the contrary in the literature. 

The computation of the retarded glueball propagator from the bulk follows the same logic 
as previously. One takes the corrected background and considers linearised fluctuations of 
the dilaton or axion about this background. The fluctuation equations will be higher order 
and will therefore have more than two independent solutions. However, only two of these 
solutions will deform into the uncorrected solutions in the a' — > limit. These are the modes 
dual to the glueball operator. The other modes are stringy degrees of freedom arising in 
the a' corrected action, these are dual to other boundary operators. 

The boundary correlator is then calculated as usual by considering the on shell corrected 
bulk action as a function of the boundary field values. As we tend towards the boundary, 
u — > 0, the leading u dependence of the field must be the same as in the uncorrected case, 
otherwise the a! corrections would not be subleading and the perturbation expansion would 
not be reliable. Thus we have from (22) 

<f> u (u) = A{uj) (1 + •••) + B{lo) (u 2 + •••) as u^O. (25) 

The coefficients A{oj) and B(u) are fixed by the incoming wave boundary condition at the 
horizon (21). The fact that the u dependence near the boundary is as in the uncorrected 
case, although the u dependence of (25) is of course different, implies that the Gr(oj) 
correlator has the same structure as previously in (20) 

Gr{u) ~ — — — ~ ~TT~\ ' ( 26 > 

<t>u<t>l \u=u B ->0 A (") 

This expression should be multiplied by functions with u dependence, and the evaluation 
as u — ► will require regularisation. By dividing by w </>* we have explicitly imposed 
the additional normalisation 4>(ub) = 1 at the boundary cutoff. The expectation that the 
result reduces to a ^<9«^>* term is realised explicitly in [29] who consider a partial set of a' 
corrections to the Yang-Mills shear viscosity. 



16 



The expression (26) was a basic connection between bulk quasinormal modes and poles 
in the dual propagators observed in [11]. A quasinormal mode in asymptotically AdS space 
satisfies incoming boundary conditions at the horizon and is normalisable as u — ► 0. This 
will occur at frequencies where A(u) = 0, which hence correspond to poles in Gr(uj). The 
corrected fluctuation equation will have quasinormal modes and hence dual poles that are 
simply deformations of the strict strong coupling poles. However, the nontrivial question 
which we will now address is whether these poles can furthermore become branch points 
and whether extra branch points can appear. We are able to give a definite negative answer 
in some cases, and expect this result to hold in general. 

Branch cuts in the u plane may appear in A{u) and B{uj) as defined in (26). The 
analytic properties of these functions are inhereted from the analytic properties of 4>oj(u) 
for small u as a function of u. Recall that u specified the boundary conditions asu^ 1. 
We shall shortly adapt some arguments from [30] to prove some analyticity of (j)^ (u) in oj 

-2-7T < 3mw =>• analytic. (27) 

Thus we gain a strip in the lower half plane in which <f>u,(u) is analytic and hence Gr(lu) 
does not have branch points. However, it was found in [12] that as p becomes large, the 
imaginary part of the location of the lowest quasinormal pole moves towards the real axis. 
In particular, there exist infinitely many values of p for which the lowest quasinormal pole is 
within the range (27). The location of these poles will only be altered by amounts of order 
a' and so will remain in this strip. However, there can be no branch points in this region 
and therefore these poles cannot have become branch points and no new branch points can 
have appeared in this region. This is the desired result. There does not seem to be a reason 
why these larger p poles should be special at strong coupling, so we expect that the absence 
of branch cuts holds in general. 

We first prove analyticity of </> u (u) in the region (27) for the uncorrected potential. 
It is convenient to work in Regge- Wheeler coordinates, so that equation (22) assumes a 
Schrodinger form. Set 

dr * = (28) 

which implies the range < r* < oo. Now let 

^ = tr 3 / 4 0. (29) 

The bulk equation (22) becomes 

-^ + V(n(u))iP = LU 2 iP, (30) 
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where 



7T 



(3u 2 + 1 



i2l 



(31) 



V(n(u))= P 2 (l-u 2 ) + — 
The only property we need of this potential is the exponential decay towards the horizon 

V{r*) ~ e' Anr * as r* ->■ oo . (32) 



The incoming boundary condition at the horizon (21) becomes 



tpuj ~ e* wr * as — > oo . 



(33) 



Now rewrite the Schrodinger equation (30) as an integral equation incorporating the 
initial condition (33) 



1 f°° 

ipu>(r*) = e lur * / cfasinu;(r* — x)V(x)ip ul (x) . 

u Jr* 

The integral equation may be solved using a Born series 

oo 

where 



(34) 



(35) 



n=0 



1 Z* 00 

-— / dxsmu(r* — x)V (x)tp^ l ~ 1 \x) 

w Jr. 



Using the inequality 



|sinu;(r* - x)\ < e -\**«\(r*-*) for x > r* , 



(36) 



(37) 



it is straightforwardly shown that 

»=o' LM ^* 

Therefore the Born series is absolutely convergent if 

dx\V(x)\e x{l3m ^~ 3miv) < oo. 



(38) 



/ 



(39) 



Given the exponential decay of the potential (32), it follows that the integral (39) is finite 
in the region (27). The consequent absolute convergence of the Born series establishes 
continuity of ip^ in this region. Analyticity follows from performing an entirely analogous 
computation for the derivative of ip^ with respect to lo. Thus we prove the analyticity claim 
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above for the uncorrected case. Whilst the lower limit of integration, r*, is small, we may 
keep it finite and therefore do not need to worry about divergences as — > 0. 

To extend this result to include a' corrections we need to perturb the solution and the 
equation. In terms of the solution we have 

ip -> V> + ei// , (40) 

where e is a small parameter controlling the a' corrections. The zeroth order solution tp 
satisfies (30) whilst the perturbation satisfies 

-^ + y(r*(«))V>' + 0^ = wV. (41) 

where O is a higher order differential operator which contains the effects of the a' correc- 
tions. This equation is the same as that satisfied by the zeroth order solution but with an 
extra inhomogeneous term. We do not need to know the precise form of O, but note that 
consistency of the a' expansion requires 

\Vi/j'\,\0^\ < -\Vil>\ as r*^oo. (42) 

The integral equation satisfied by the perturbation is 

1 f°° 

^(r,) = e lujr * - - / dxsmoj{n - x) [V{x)^{x) + Cty w (z)] . (43) 

We may again solve this equation using a Born series. All the terms are as previously in 
(36) except for the first term which is now 

i r°° 

V4 0) (r*) = e^ r * - - / dx sin u;(r* - x)Oij> u {x) . (44) 
The zeroth order solution, ip^^x), satisfies the inequality 

\Mx)\<te- xJmu , (45) 

where k is some constant. The bound (45) follows from the asymptotic behaviour of i/) w (x) 
(33) and the absence of divergences at finite values of x and as x — > 0. Using this bound 
together with the inequalities (42) it follows that 



k£ n) (r*) < e-^l^^l ( [1 + 7] exp / dx\V{x)\e x ^ Jm ^-° m ^ -M . 

n=0 L I I r* J / 



(46) 



Thus if (39) holds, the Born series is absolutely convergent. Along with an analogous result 
for the derivative of ip'^ with respect to u, the convergence again implies analyticity of ip' w 
in the strip (27). This completes the proof of our claim including a 1 corrections. 
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5 Discussion 



We have seen that frequency space glueball propagators in large N thermal N = 4 SYM 
theory on R 1,3 have different analytic structures at weak and strong 't Hooft coupling. 
In position space this difference translated into qualitatively different time decays of field 
fluctuations about equilibrium. In the perturbative regime the correlators are known to 
exhibit power law falloff on intermediate time scales followed by a slow exponential decay 
at late times. We contrasted this with the rapid pure exponential decays seen at strong 
coupling. 

Based on our free field computations and the generally understood aspects of weakly in- 
teracting plasmas, we expect that the perturbative real time response functions will exhibit 
branch cut non-analyticities. These differ qualitatively from the strong coupling quasinor- 
mal poles, which we further argued survive away from the strict infinite coupling limit. 
This picture suggests that these propagators must depend non-analytically on the 't Hooft 
coupling. The non-analyticity could either be due to nonperturbative corrections to both 
weak and strong coupling limits or else to a phase transition at an intermediate coupling. 
As we noted, however, further work is needed to rule out the possibility that the branch 
cuts seen in perturbation theory actually split into a series of poles in the large N theory, 
thus allowing a continuous interpolation to strong coupling. 

One observation that may favour the phase transition possiblity is that to each order 
in the large N expansion, the perturbative 't Hooft expansion for zero temperature N = 4 
SYM theory is thought to be regular at A = 0, with a finite radius of convergence. This is 
due to the slow growth of the number of planar diagrams at each order in A [31]. Including 
a finite temperature should not modify this regularity beyond introducing a dependence on 
y/X. The convergence of the 't Hooft expansion does not rule out nonperturbative terms 
in A, but it does mean that they are not necessary, unlike in standard Yang-Mills theory. 
If the non-analyticity is not due to nonperturbative corrections, then our results predict a 
phase transition in finite temperature Yang-Mills theory on M 1,3 . 

Regarding corrections to the strong coupling results, we have argued that perturbative 
a' corrections do not lead to a broadening or merging of the poles into branch cuts. In 
principle there could be nonperturbative a' corrections to the strong coupling results, but 
there are no obvious worldsheet instantons in the bulk. Our strong coupling considerations 
have all been in the strict large N limit. It is possible that the quasinormal poles are 
broadened into branch cuts by 1/N effects and it is expected that 1/N corrections will 
introduce additional branch cuts in real time correlators at all couplings [24]. 
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A natural candidate for a dual phase transition in the bulk is given by the Horowitz- 
Polchinski correspondence [32, 33]. As the AdS background becomes increasingly curved, 
the horizon size eventually becomes of the order of the string scale. The Horowitz-Polchinski 
correspondence suggests that in this regime the appropriate description of the system is as 
an ensemble of highly excited string states. It would be interesting to ascertain whether 
there is a connection between the qualitative change in behaviour we find in field theory 
from strong to weak coupling and the putative change in the bulk from a black hole to 
string states. 

These possibilities have been entertained before. A phase transition at a critical 't Hooft 
coupling in finite temperature N = 4 SYM theory was previously argued for in [34], mainly 
on the basis of analytic properties of the free energy as a function of coupling. It was then 
suggested that this phase transition was dual to the Horowitz-Polchinski point in [35] 2 . A 
key question for placing the arguments about phase transitions on a firmer footing is the 
identification of an appropriate order parameter. One would also want to establish whether 
the phase transition existed only in dynamical response functions or whether the static 
equilibrium theory was also affected. 

Finally, in a different direction, various proposals have been made for isolating features 
of the bulk black hole singularity in dual field theory correlators [36, 37]. These works have 
focused on infinitely massive bulk probes in order to employ a semiclassical approximation. 
The correlators we have computed in field theory are dual to massless minimally coupled 
bulk scalars. Nonetheless, it would be interesting to see if signatures of a singularity, possibly 
resolved, can be seen in our expressions. 

Acknowledgements: During this work we have had useful discussions and correspon- 
dence with Nick Dorey, Michael Green, Shiraz Minwalla, Asad Naqvi, Toby Wiseman and 
Larry Yaffe. We are especially grateful to Ofer Aharony, Dan Boyanovsky and Rob Pisarski 
for critical comments and extensive discussions. SAH is supported by a research fellowship 
from Clare College, Cambridge. SPK acknowledges support from a PPARC Rolling Grant. 

2 This paper further argued that the Hawking-Page transition would not continue to weak coupling; but 
this now appears not to be the case [2, 3]. However, that argument seems to be independent of the arguments 
in the same paper and in [34] for a phase transition as a function of the 't Hooft coupling. 
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A Effect of thermal mass and width 



We discuss how the inclusion of the thermal mass my ~ \f\T and thermal width V for 
the intermediate gluon states in the correlator (6) changes the analytic structure of the 
frequency space propagator. Including a mass changes the frequencies in the real time 
propagators (4) to E = uj p = p 2 + m^. Restricting attention to the homogeneous case 
(the inhomogeneous case is similar) the retarded glueball propagator becomes 

(G°(t)g°(0)) R = -4^0(t) p' 2 d / 2p ' 2 + f T)2 coth (^) sin(2uyi) , (47) 

where Ccy = \jp' 2 + m^. Upon Fourier transforming to frequency space this function is seen 
to have branch cuts. The easiest way to see this is to note that upon Fourier transforming 
&(t) sin(2uyi) using an ie prescription we get 

j2 />oo (^Jl i ™2^2 ,, , s r 1 1 

• (48) 



G° R (u) = 2— / p' 2 dp' y v 2 T> coth 



w — 2cj p / + cj + 2oj p / + ie 



p' 

The presence of discontinuities can now be inferred from the imaginary parts of the energy 
denominators in the above expression by using the representation 

— nr- - = p — A — + mt<5(w ± 2bv) . (49) 

lo ± 2uj p i + ie L0 ± 2uJ p r 
It is now a straightforward excercise to see that (48) has a square root branch cut discon- 
tinuity between u> = Itxit and oj = —2rriT- Note that this cut on the real axis would also 
appear for a massive particle at zero temperature. 

Finally, we turn to the inclusion of the thermal width in the propagators. This is a subtle 
issue since it involves unstable intermediate states. We can provide a heuristic argument 
for its effect by simply looking at the form of the propagators for a massive but unstable 
resonance at zero temperature. The propagator for such a resonance, at zero temperature, 
is (pq —p 2 — m 2 + where V is the width of the resonance. Repeating the above excercise 
in the real time domain and computing the frequency space retarded propagator, it is clear 
that the branch cut along the real axis gets shifted off the real axis. Whether this conclusion 
is correct depends on a more detailed and precise field theoretical analysis of the analytic 
structure of self-energy corrections to the internal lines. 



B Thermal glueball propagators on R x S' s 

In this appendix we indicate how our field theory computations may be generalised to 
thermal Yang-Mills theory on 1 x S 3 , The glueball time dependence in this context is 
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interesting as it may be a formalism in which one can study nonequilibrium processes such 
as black hole evaporation. 



B.l Expansion in harmonics 

Free Yang-Mills theory on K x 5 3 has a nontrivial phase structure because the finite spatial 
volume restricts the partition function to a sum over colour singlets [2, 3]. As emphasised 
in [3] this may be seen in the fact that the nondynamical mode 

is strongly coupled even in the g 2 N — > limit and so needs to be treated nonperturbatively. 
The different phases are characterised by the distribution of eigenvalues of a [3] . 

Working on S* 3 x R, it is useful to decompose the gauge potential into transverse vector 
spherical harmonics on S 3 

A?(t,e) = J2A am Wr(0), (51) 

m 

In this notation, m denotes all the relevant quantum numbers for the harmonic decompo- 
sition. For scalar harmonics these are {j m ,m m ,n m } and there is an extra e m = ±1 for the 
vector harmonics [4, 38]. The quadratic part of the Yang-Mills action then becomes 

C 2 = ~A a ™ [D 2 t + E 2 m ] A am . (52) 

As previously, in free field we may work directly with the physical transverse degrees of 

freedom and ignore the other modes and ghosts. In expression (52) 

ft Aam 

D t A am = + f abc a b A cm , (53) 

and the energy levels are 

E m = j m + 1 • (54) 

We are working in units in which the radius of the S* 3 is one. 

The glueball operator (2) may also be expanded into harmonics. We restrict ourselves 
here to the TrF 2 operator. The quadratic term gives 

Q 2 = 2 [-D t A am D t A an D mnp + e m e n E m E n A am A an D mnp ] S p . (55) 

This expression should also contain terms dependent on the time component Aq of the 
gauge field, but these do not contribute to the free theory propagators. The expression (55) 
is obtained using completeness of the harmonic decomposition and the triple overlap [38] 

D mn P= f d Q V im v n S p^ ^ 
JS3 
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where S and V are the scalar and vector harmonics, respectively. 

The thermal propagator following from the action (52) is now, for t > 0, 

jp—fat 

(A am (0)A bn (t)) = S mn —r- [(n(E m ) + \)e~ iE ^ + n{E m )e iE ^] , (57) 

-'Arc 

where we have suppressed the colour indices and used the abbreviation 

{far = f abc a b (5g) 

The Bose-Einstein distribution is now 

n(E) = alw \ , . (59) 

v I e /3(E-ifa) _ l ^ ' 

We can see in (59) that the a condensate plays the role of a chemical potential. This 
corresponds with the fact that a is a Lagrange multiplier for the conservation of colour 
charge. 

B.2 The free theory two point function 

We find that for t > the two point function of the glueball (55) is given by 



Qmn(t) 

1+ I- 



(60) 



where the first sum runs over 

1+ = {m,n\ \j m ~ 3n I < jp < jm + 3n - 2} , (61) 



and the second over 



I- = {m, n | \j m - jn\+2<j p < jm + jn} ■ (62) 



The coefficients are 



<•-„„ [{E m + E n + E v f - 1] [{E m + E n - E v f - 1] , 



^mnp 



c mnp = [(E m + E P - E n f - 1] [{E n + E p - E m f - 1] , (63) 

and the time dependence is contained in 

Q+ n (t) = (n(i? m ) + l)(n(i?0 + l)e-^ + ^ )i +n(^ m )n(S n ) e ^ + ^) i , 
Q~ n (t) = (n(E m ) + l)n(E n )e-^- E ^ + n(E m )(n(E n ) + l)e^- E ^ . (64) 
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In deriving these results one needs to use identities for 3j-symbols that may be found for 
example in the appendices of [4]. 

The trace in (60) is over the colour indices of fa as given in (58). This may be expressed 
in terms of the eigenvalues Aj of a considered in [3] as follows: For some function F of fa 
we have 

N 

(65) 



TrF(fa)=J2F(X i -X j ). 

In practice, we need to evaluate this sum in the two possible vacua identified by [3]. The 
first is the maximally clumped vacuum in which all of the TV eigenvalues are equal 



Aj = const. 



TrF(fa) = N 2 F(0) . 



(66) 



This is the correct vacuum at high temperatures. The other possibility is a uniformly 
distributed vacuum 



2iri , , %-~ I ( 2irn \ / n \ 

— =► TrF(fa) = 2N^F( — )(l--) + NF(0), 

n=l ^ ' 



A, = 



(67) 



which is the correct vacuum at low temperatures. In this last expression we assumed that 
F is an even function which is indeed true for (60). 

The two point function (60) has no mixing between the different <S 3 harmonics, as we 
expect for a free theory. The simplest case to consider is the propagation of the homogeneous 
mode p = q = 0. In this case from (62) we have E m = E n in (60) and only the first sum 
contributes to the time dependence. For the \ = const vacuum the sums in (60) for the 
homogeneous mode admit a compact expression in terms of theta functions 

-N 2 \ .d d 3 l d ^[{t,e-^ 2 ) 



&<0°(O)S°(t)> 



2vr 2 

—N 2 
~2^ 



A d d 3 

—4 1 

dt ^ dt 3 



dt 2 



dt 4 



dpti^e-P/ 2 ) ' 

'-&[(t,e-^ 2 ) ^(t,e^ 2 )' 
^i(t,e-/3/ 2 ) ~ «?i(t,e^/ 2 ) 



(68) 



These expressions are true for t > 0. We have dropped various S(t) contributions to the 
real part. We have not attempted to generalise (68) to the inhomogeneous modes. 

One significant feature of the result given in expressions (60) to (64) is that it is periodic 
in time. The periodicity of (60) under t — > t + 2ir is a special case of the quasi-periodicity 
that one always finds for field theories on a compact spatial manifold. The exact periodicity 
in this case arises because in our free theory on 5 3 x K the energy eigenvalues are evenly 
spaced. The quasi-periodicity is expected to be dual to large gravitational fluctuations in the 
bulk [39] such as the formation and evporation of black holes. However, this correspondence 
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has only been made precise when averaged over long times [40, 41]. The requirement of long 
time averaging is due to the thermal nature of the semiclassical approximation in gravity. 

At increasingly high temperatures, f3 — ► 0, the nontrivial time dependence is restricted 
to times t < /3, much shorter than the full period. This limit, which can be viewed as letting 
S* 3 — > M 3 , is simpler to study analytically was the subject of the prvious three sections. Note 
however that the nontrivial phase structure uncovered in [2, 3] occurs at finite temperatures. 
In particular, addressing black hole evaporation in this formalism will require working at 
finite temperature. 

At low temperatures, — > oo, we need to use the vacuum with evenly spaced eigenvalues 
[3]. Thus we use the formula (67) to evaluate the traces. In particular, evaluating the leading 
finite temperature correction will require the result 

Trcos/3/a = ^J. (69) 



sin ^ 



The resulting two-point correlator is then 

T2 



WW)) = - -^57 [cost + 2 cos 3 t] - -p-^-e-We-™ + • • • . (70) 
vr 2 sm° t 7r 2 sm 2 ^ 

The retarded propagator is then 

(g(0)g(t)) R = [cos t + 2 cos 3 1] + 2e~ 2/3 %- sin 4t + • • • . (71) 

vr 2 sin ij vr 2 sm 2 ^ 

The leading finite temperature correction is exponentially suppressed. 
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